Longitudinal spin waves in a dilute Bose gas 
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We present a kinetic theory for a dilute noncondensed Bose gas of two-level atoms that predicts 
the transient spin segregation observed in a recent experiment. The underlying mechanism driving 
spin currents in the gas is due to a mean field effect arising from the quantum interference between 
the direct and exchange scattering of atoms in different spin states. We numerically solve the spin 
Boltzmann equation, using a one dimensional model, and find excellent agreement with experimental 
data. 
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A recent experiment at JILA has displayed remarkable 
effects of spin density fractionation in a trapped, ultra- 
cold gas of Rb atoms with no Bose-Einstein condensate 
present G[. Under conditions which wc summarize in 
more detail below, sudden preparation of all atoms in 
a coherent superposition of two spin states generates a 
spin wave resulting in the observed spatial separation of 
the two components. This occurs even though both the 
mean field and differential Zeeman energy differences are 
almost a thousand times smaller than the thermal energy 
kgT. In this paper, we show that this astonishing depar- 
ture from equilibrium results from quantum interference 
between direct and exchange scattering of atoms in the 
two spin states. A first-principles kinetic theory with 
no fitting parameters gives excellent agreement with the 
experimental data, and suggests possibilities for quan- 
titative studies of quantum coherence in noncondensed 
Bose-Einstein gases. 

In the experiment of Lewandowski et ai, a few mil- 
lion atoms of 87 Rb are confined in a magnetic cigar- 
shaped trap (uj z /2tt = 7 Hz, uj r j2-K = 230 Hz). By 
applying microwave and radio-frequency radiation, all 
atoms in the gas can be uniformly prepared in an ar- 
bitrary superposition of the \F = l,Mp = —1) = |1) 
and |2, 1) = |2) hyperfme states of the ground configu- 
ration. The frequency splitting, A = u>\ — 0J2, between 
the two states depends on the position, r of the atom 
in the trap: A(r) = Abr(i") + Amf(f) — Whf, where Whf 
is an overall uniform frequency splitting (whf/27r ~ 6.8 
GHz). The first term Abr(i*) is due to the differential 
Zeeman effect, predicted by the Breit-Rabi formula, for 
atoms in a nonuniform magnetic field. The second term 
is due to the mean field frequency shift proportional to 
the density of the gas, which has a Gaussian profile in the 
harmonic trap. By applying two 7r/2 pulses separated by 
a variable delay time, this local frequency splitting is ex- 
tracted from the Ramsey interference fringes measured at 
different positions along the axial direction of the trap. 

Lewandowski et al. describe experiments in which the 
second 7r/2 pulse is omitted and the time evolution of the 
density of either state is observed after the initial 7r/2 
pulse. The following spectacular behavior is observed: 



the densities of the two states segregate along the ax- 
ial direction of the trap and then relax to a completely 
overlapping stationary state after approximately 200 ms. 
This fascinating behavior is found to depend crucially on 
two different parameters: the density of the gas n and 
the nonuniformity of the local frequency splitting A(r). 
No segregation is observed when the density is lowered 
below a critical value, or A(r) is made approximately uni- 
form by adjusting the bias magnetic field. On the other 
hand, the segregation effect becomes more dramatic as 
the density and the inhomogeneity of the splitting are 
increased. 

In this letter, we show that the transient spin segrega- 
tion is actually an overdamped spin wave arising solely 
from the mean field of the gas, even if the interaction is 
spin independent. This effect is well known from earlier 
work done on spin polarized hydrogen gases J|, |[ ||, |), 
and is initiated by the spatially varying local frequency 
splitting A(r). When A(r) is uniform, the spins of the 
atoms throughout the gas precess in exactly the same 
fashion so that every forward scattering event can be un- 
derstood in terms of two identical atoms colliding, giving 
rise to the well-known factor of 2 from the direct and ex- 
change terms in the Hartree-Fock mean-field theory of a 
Bose gas. However, when A(r) depends on position, two 
colliding atoms will have acquired different spin states de- 
pending on their history in the trap. In this case, when 
the direct and exchange forward scattering events are 
added, an additional mean field term appears propor- 
tional to the local spin S(r, t) of the gas, which accounts 
for the constructive and destructive interference between 
the two scattering paths. It is this term that gives rise 
to the spin wave. This mean field effect occurs when the 
transverse spin (the internal coherence) is long lived com- 
pared to thermal relaxation, emphasizing that the gas of 
two-level atoms in the JILA experiment is quite differ- 
ent from an incoherent binary mixture. A solution of the 
collisionless Boltzmann equation for the spin is already 
sufficient to predict spin waves. However, because the 
JILA experiment is in an intermediate regime approach- 
ing the hydrodynamic region, the spin current is strongly 
damped due to collisions. 
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The Hamiltonian describing a single, trapped, two- 
level atom of mass m is: 



The collision integral in Eq. (||) is given by 



H = 



V 2 + t/ cxt (r) 

2m 



(1) 



The first term in ([!]) is the center of mass Hamiltonian 
containing the kinetic energy and the external parabolic 
trap £/ C xt(r) = mu 2 [a 2 (x 2 + y 2 ) + z 2 ] /2, where a = 
Lu r /to z . This part of the Hamiltonian is uncoupled from 
the internal, pseudo-spin, degree of freedom, which is 
governed by the second term: O(r) • f = O x (r)f a + 
Qy(r)f y + Q, z (t)t z , where Tj is a Pauli matrix. In the ab- 
sence of an external coupling field, Q z — Abr (we make 
the rotating wave approximation to eliminate the hyper- 
fine splitting Whf )■ We model binary interactions between 
particles by a delta pseudo potential describing elastic, 
spin preserving collisions, the strength of which depends 
on the hyperfine states T4,(r,r') = gijS(r — r'), where 
gij = Ann a,ij/m, with a%j being the scattering length 
for collisions between atoms of species i and j. For 87 Rb, 
we take an = 100. 9ao, a%2 = 98.2ao, 022 = 95.6ao, where 
ao is the Bohr radius G[. 

Several groups have previously worked out the fun- 
damental kinetic theory of a noncondensed dilute Bose 
gas with internal degrees of freedom, to describe spin 
waves in spin-polarized atomic hydrogen || |^, [7], ^, |^, 
|lO| , O , p^f . Using a semiclassical approximation to de- 
scribe atomic motion in terms of a phase space distri- 
bution function, we obtain coupled Boltzmann equations 
for the distribution functions of atomic density, /(r, p, t), 
and spin density, cr(r,p,i): 



ot m 



VU, 



■v p /--vn r , 



• V„(T„ = — 



dt 



—+~Va-VU„-V p a--Vn n -V p f~n n xa = — 



coll 

(2) 



coll 

(3) 



Eq. (H) has an implicit sum over the repeated index a 
The total density and spin density are obtained from the 
distribution functions as n(r,t) = ni(r,t) + n 2 (r,t) = 
J dp/(r,p,i)/(27r?i) 3 and S(r,t) = J dpa{r,p,t)/(2nh) 3 
respectively. Here the longitudinal component of the spin 
represents the relative density S z — ri\ — rii and the 
transverse components S x and S y describe the real and 
imaginary parts of the internal coherence. The center of 
mass effective potential is U n (r,t) = U cxt (r) + gntii + 
522^2 + 312(^1 + "2)/2. The modified coupling field in- 
cluding mean-field effects is 



(4) 



where n' n (r, t) = {£l x (r), tt y (r), ft z (r) + AmfO, t)}, and 
A M f(i", t) = 2 [gnni + gi 2 n 2 - (.922"2 + 912^1)] /Ti. (5) 



da 
dt 



dp 2 



dp 3 



^12 

h J (2irh) 3 J (2-Khy 
5{e p + e p2 - e P3 - e P4 )5(p - 
x{3/(P3)^(P4) + CT(p 3 )/(p 4 
-/(p)ff(pa) - 3c?(p)/(p 2 )}, 



dp 4 

P2 P3 P4) 

(6) 

where e p = p 2 /2m. Here we neglect a principal value 
contribution, which gives a second-order correction to 
the free streaming evolution, and we take all scattering 
lengths <2y to be equal - a reasonable approximation for 
87 Rb. This approximation results in the conservation of 
spin density during collisions, i.e. J dpda/dt\ co u = 0. 
When the small differences in scattering lengths are ac- 
counted for, the transverse spin decays slowly. For 87 Rb, 
this contribution to the "T2" lifetime is of the order of 
10s@. 

Immediate insight can be gained if we solve for the 
time evolution of the spin density by integrating (^) over 
momentum 



dS_ 
~dt 



1 



V • J = 0' x S. 



(7) 



Since 5x5 = 0, the second term g^S/h in (||) has no 
direct affect on the spin density. We show below that 
this term instead sets up a spin current J in the gas 
that strongly affects 5. Also, an interesting paradox has 
emerged concerning the factor of 2 in the mean-field fre- 
quency shift. For an incoherent binary mixture of atoms 
in either of the states |1) or |2) (i.e. a x — a y = 0), 
it is straightforward to show from (Q) and (0) that the 
difference in chemicalpotcntials due to interactions is 
Ml - M2 = [2(5n?ii - 522«- 2 ) + 512 {ri2 - ni)}. There is 
a factor of 1 instead of 2 in front of g\ 2 since the two 
atoms are distinguishable. However, when the atoms are 
manipulated in a coherent fashion, from (j7j) we see that 
the precession of the transverse spin is given by Amf, 
which is the mean-field frequency shift measured by the 
Ramsey interference technique. 

The hydrodynamic equation for the spin current J is 
obtained by assuming the following simple form for the 
spin distribution function a(r,p,t) = /o(r, p)[M(r, t) + 
p ■ v(r,£)/fc B T], where M(r,i) = S(r,t)/n (r) is the 
renormalized spin density. Here, and in the rest of 
this letter, we assume the that the total phase space 
density is stationary |l6| / (r, p) = 770 exp{ — [p 2 /2m + 
U cxt {r))/k B T}. The constant 770 = a 2 N {Twj z / k B Tf is 
determined by the requirement J dp J dr fo / (2irh) 3 = N, 
where N is the total number of atoms. The equilib- 
rium total density is no(r) = J dpfo(r, p) / (2irh) 3 . Using 
this ansatz, the equation of motion for the spin current 
J = J dppa j \2irh) 3 is 

^-0„X J + k B TVS + SVU+ ^n Vn n = — . (8) 
dt 2 td 

Here, the diffusion relaxation time is [O, [l^] rc(r) = 
[(32a 2 2 n (r)/3)^Trk B T/m]-\ 
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The spin segregation dynamics is described by the lon- 
gitudinal spin S z . Taking the z component of (fjj) and (||) 

nves 



dS z 
dt 



1 



3J, 

dt 



V J 2 



n t d 



(9) 



(10) 



where F z 



-k B TVS z - S Z VU - hn a Vn nz /2. In the 



absence of a coupling field fl x — fl y — 0, the evolution of 
S z is entirely due to the spin current J z and thus the total 
S z is conserved J dzS z = 0. The current 3 Z is driven 



by the first two terms on the right-hand side of Eq. (10). 
The first term is due to the mechanical force F arising 
from the spatial gradient of the local energy splitting. 
The second term represents the spin current driven by 
the dynamics of the S x and S y components through the 
term gi^S/h. The third term represents the diffusion 
transport process that gives rise to the damping of the 
current. As already mentioned in this paper and also 
discussed in Ref. M, the magnitude of the mechanical 
force F is negligibly small for driving the spin current J z 
in the experimental situation. To highlight the effect of 



the term g^S/h, we take the time derivative of Eq. (|10[), 
work to first order in gi2, and neglect F 2 , the relaxation 
term, and terms of second order in J. This gives 



d 2 J z 

dt 2 



9i2k B T , ^ ^ Sx giiknT 2 



-(SxVS) z = --- 



(11) 



where S± and <j) are the amplitude and phase angle 
of the transverse spin component S x + iS y = S±e 1 ^. 
Eq. ( pd| ) explicitly shows that the spatial gradient of 
the phase angle induces the spin current J z . In a 
short period of time right after the ir/2 pulse, it is 
reasonable to assume that the transverse spin compo- 
nents are undergoing rotation with the local Larmor fre- 
quency 4>(r, t) « n nz (r)t. With this simple approxima- 
tion, the induced spin current for short times is given 
by 3 Z ~ -g 12 k B TSl(r)Vn nz (r)£/6h. Taking fl nz = 
Abr + Amf ~ — muj 2 iS z 2 /2h |l7f| and the initial con- 
dition 5_i_(r) = n o (0)e _l/cxt / fcBT , we find that the initial 
evolution of S z after the it/2 pulse is given by. 



S z (r,t) 



gi 2 n§(r)fc B r 2 



1 



2mu 2 



t 4 



(12) 



The above form ula predicts that S z has nodes at z = 
±y // k B T/2mcu 2 , which we have verified in the numerical 



calculation below. Eq. (12) identifies the characteristic 
timescale i sp i n = [2Ah 2 / gi2n (0)k B Tuj 2 iS ] 1 ^' 1 needed for 
the system to build up the S z component. For the JILA 
parameters, t spm w 30 ms, which is consistent with the 
delay time seen in experiments before the segregation 
begins. 

To complement our simplified analysis, we also numer- 
ically solve the Boltzmann equation (|^). Motivated by 



the observation that spin segregation occurs only along 
the axial direction [jjj, we construct a one dimensional 
model of the system by making the ansatz a(r, p, t) = 
a(z,p,t)ho(x,y,p x ,p y ) and then averaging over x and y. 
Here we take the static profile in the radial direction to be 
of Gaussian form h = cxp{ — [(p 2 + p 2 )/2m + mwj(i 2 + 
y 2 )/2]/k B T}. We substitute this ansatz into (||) and in- 
tegrate over the radial phase space variables, which gives 
the following one-dimensional model Boltzmann equa- 
tion 



da 
dt 



p da 

m dz 



dU cxt da 
dz dp 



Q n x a = 



da 
dt 



(13) 



ID 



Here we have made the approximation U n — C/ cx t 
and we have dropped the fourth term in (||) coupling 
the center of mass and spin. The collision integral 
in one dimension involves a phase space average in 
the radial direction da/dt\ 1D = J xy da / dt\ coll / J xy h , 

where we have introduced the notation f ■ ■ ■ = 

Jxy 

J dx J dy J dp x J dp y ■ ■ ■ /{2ttTl) 2 . The radial averaging 
introduces a scaling factor in the mean field terms, so 
that gij — > g[j = gij / {2\ 2 m ), where AdB is the thermal de 
Broglie wavelength, g'^ has the correct units of energy 
times distance required in our one dimensional model. 

Although the direct numerical simulation using the full 
expression for the one dimensional collision integral de- 
rived from Eq. (||) is technically feasible, we introduce a 
simple model for the relaxation 



da 
~di 



1 



ID 



(z) 



[a(z,p,t)-M(z,t)f (z,p)}, (14) 



where t c \{z) = [l6a 2 2 n o( z ) V '^k B T / m] 1 is the radially 
averaged mean collision time, fo(z,p) = /o(r,p)//io, 
and M(z.t) = S{z,t)/n (z). Eq. (Q) contains the 
essential properties of collisions: (i) it vanishes when 
the distribution function has the local equilibrium form 
a(r, p, t) oc e~ p / 2mi:BT , (ii) it conserves the spin density. 
We note that the form (|l4j) does not require the knowl- 
edge of the long-time equilibrium solution for S(r,t). 

We solved the one-dimensional spin kinetic equation 
( |l3| ) numerically using a finite difference scheme. For 
the initial state of the spin, we take a y (z,p,t — 0) = 
fo(z,p),a x (z,p,t = 0) = a z (z,p,t = 0) = 0, correspond- 
ing to the state immediately following the first tt/2 pulse. 
Figure 1 shows the time sequence of the density of the 
|1) state ni(z,t) = [n (z) + S z (z,t)]/2, corresponding to 
Fig. 3(c) columns (v)-(vii) of Ref. 0. This shows that the 
spin segregation vanishes when the density is lowered to 
n~2x 10 12 cm" 3 . The red curve is the raw JILA data 
and the blue dashed line is the theory. The agreement 
is striking. We also compare our numerical results to 
Fig. 3(c) columns (i)-(iv) of Ref. [Q, shown in Figure 2. 
In this case the density is held fixed but the curvature of 
the energy splitting is varied in each sequence. Note that 
when Abr(z) is chosen to approximately cancel Amf in 
the second column, the amplitude of the spin wave is es- 
sentially zero. For the larger values of Wdiff considered in 
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FIG. 1: Time sequence of ni(z,t). The red line is the un- 
smoothed JILA data and the blue line is the numerical solu- 
tion to (|l3|). We have taken a temperature of T = 850 nK Jl3| |, 
and the peak total density n listed under each column is in 
units of 10 13 cm -3 . We approximate Abr = ma> BR z 2 /2?i, and 
take ujbb./u> z = 0.09 for each column in our calculation. This 
corresponds to a value Wdiff/^z = 0.1 for the first column. 
The axial position z is in the range ±600 /J,m. 



FIG. 3: Modulation of Ramsey fringes. The top row corre- 
sponds to the third column of Figure 1, and in the bottom 
row we have set Abr = (which corresponds to sitting at 
the "magic spot" bias field of 0.323 mT). The first column is 
taken at z = 7 /im and the second column at z = 350 /im. The 
red and blue lines compare the dynamics with and without 
collisional relaxation, respectively. We have taken a detuning 
of 5/luz = 5. 



t=10 




v = 0.09i 
diff 



0.00 



0.16 



0.20 



FIG. 2: Time sequence of m(z,t). We have taken a temper- 
ature of T = 850 nK 113] , and the peak total density n for 
each column is 1.82 x Hr 3 
units of Hz. 



cm . Here, va\a = cJdiff/27r is in 



columns three and four, the agreement is only qualita- 
tive. 

We finally investigate the effect that the spin wave 
and relaxation have on the Ramsey fringes. Ramsey 
fringes can be obtained from the results of our calcu- 
lation by simply rotating the Bloch vector S(z,t) at 
each time by 90° about the local oscillator vector SI = 



{cos(<5f), sin(Jf ), 0}, where 5 is the detuning between cou- 
pling field and the hyperfine splitting. In the top row of 
Figure 3 we show the Ramsey fringes corresponding to 
the third column of Figure 1, taken at the center (left) 
and the edge (right) of the cloud. We compare the Ram- 
sey fringes with (red) and without (blue) relaxation (i.e. 
da/dt\ 1D = 0). We find that the Ramsey fringes taken at 
different positions across the cloud are all modulated at 
the period of the spin wave, and that the fringe visibility 
decays when the effect of collisions is included through 
d<r/dt\ 1D . At long times the system relaxes to the state 

5 = 0, that is, the gas evolves to a completely overlap- 
ping binary mixture, which is reflected by the vanishing 
of the fringe visibility. In the bottom row, the differential 
Zeeman splitting is set to zero Abr = 0, which reduces 
the curvature of the local energy splitting. This has two 
main effects: the frequency of the spin wave is lowered 
and the fringes are visible for a much longer time. In a 
calculation not shown in Figure 3 for the case where the 
curvature is approximately zero, as in column two of Fig- 
ure 2, we find that the Ramsey fringes are not modulated 
and do not decay, emphasizing that collisional relaxation 
damps out spin currents, while conserving the spin. The 
trends shown in the bottom row of Figure 3 agree quali- 
tatively with experiment [l3| . 

In summary, we have presented a simplified model for 
the spin density and current that predicts longitudinal 
spin waves in the gas when the energy splitting between 
hyperfine states in nonuniform. We have also numerically 
solved a one dimensional model of the Boltzmann equa- 
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tion that supports our simplified hydrodynamic model 
and gives excellent agreement with the JILA experi- 
ment 

While this manuscript was being written, two inde- 
pendent articles 14 Eq] appeared on the lanl.arXiv.org 



e-print archive that make predictions similar to ours. 
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